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ABSTRACT 

Multi Resource Generalized Assignment Problem (MRGAP) 
is to determine an optimal allocation of tasks to agents such 
that the requirements of tasks allocated to agents do not exceed 
the available resources with ag^ts* The MKJAP is a generali- 
zation of the well-known generalized assignment problem in the 
sense that each agent possesses multiple re 60 in:*ces and ^loca- 
tion of tasks must satisfy each resource limitation. We des- 
cribe several applications of MKJAP and develop a bran^ and 
bound algorithm based on lagrangian relaxation concepts. A 
novel feature of the algorithm is that the lagrangian multi- 
pliers are confuted only once and these multipliers are used 
to compute the lower boujid for each siibproblem. The proposed 
elgorithm was coded and tested extensively on randcwaly generated 
problems of different parameters. Results of these computa- 
tions are quite encouraging. The algorithm can solve problems 
with 1000 0-1 variables and 6 resources optimally in less than 
■60 seconds and problems with 8000 variables and 6 resources 
within 1 percent of optimality in less than 5 minutes. 



CHAPTER I 


INTRODUCTION 


1.1 INTRODUCTIONS 

The formixlation of models with different mathematical 
structxares and the development of algorithms to exploit these 
structures continues to be a productive direction for research 
in mathematical programming. During the last few years, several 
similarly structured integer programming models, namely, 
weighted assignment models, have been developed for a number 
of diverse applications. In this thesis, we deal with one 
such mathematical model, which comes into the category of 
weighted assignment models. Multiple Resource Generalized 
Assignment Model, which will be referred to as MRGAP, here- 
after. This problem is a known NP-complete problem and we 
develop a branch and bound algorithm for this problem based 
on lagrangian relaxation concepts. 

In the classical assignment model [lO], we have agents 
and tasks are done by these agents* The aim is to find out 
the optimal pairings of agents and tasks. The optimality is 
determined by a criterion function like cost etc. Each task 
is to be assigned to a single agent and each agent Should be 
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given only one task* These constraints restrict the model to 
a limited domain of applications like assigning air crew to 
the flights, assigning jobs to machines or men and space to 
departments etc. However, in many practical situations the 
number of tasks to be done are more than the number of agents 
and accordingly more than one task should be given to each 
agent. This leads to a more useful model which allows for • 
more than one task to be assigned to each agent, provided 
these tasks would not require more of some resource than 
is available with the agent. This type of model, which is 
given the name Generalized Assignment Problem (GAP) ty Ross 
and Soland [ll], has many applications like assigning jobs 
to coirputers in a computer network, assigning software deve- 
lofraent tasks to programmers and scheduling the variable 
length television commercials into time slots etc. 

■ From the applications stand point, a more useful 
model would be the one which allows for more than one reso- 
urce with each agent and which allows the tasks to be comple- 
ted in two or more stages, that is, the tasks will be completed 
by the agents at different levels. These situations generally 
arise in sophisticated systems where everythiiiig (for example, 
manpower, equipment, capacity, etc.) is considered to be very 
critical and has to be taken care of while finding the best 
assignments. These mpdels are temed as Weighted Assignment 
Models (WAM) in literature. 
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In this thesis, we consider a special case of the WAm 
in which tasks are to be finished at a single level, that is, 
if a task is assigned to an agent, it will be finished comple- 
tely by that agent at a time. Each agent will have more than 
one resoiorce which will be consumed when a task is assigned 
to it. The novel features of our algorithm are that lagrangian 
relaxation concepts are used to obtain tight lower bounds and 
the lagrange multipliers are calculated only at the root node 
and are used throughout the branch and bound tree. The LIFO 
treatment of candidate problems had been employed in the 
branch and bound algorithm. Numerical investigations with the 
algorithm are presented ^n detail. 

1.2 STATEMENT OF THE PROBLEMS 

The notations that are used throughout this thesis 
have been given below for the sake of completeness. 

Set of agents, 

Set of tasks 

|j|. 

Number of resources associated with each agent, 
1, if task j is assigned to agent i 
0, otherwise 

The amount of resource k available with the 
agent i. 


I 

J 

m 

n 

P 


X. 


10 

k 
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S Cost of assigning task j to agent i, 

t j^mount that will be consijmed of resoijrce k 
when task j is assigned to agent i. 


In this thesis, we have considered a model where a 
set of tasks must be divided among a set of agents and each 
task must be completed by only one agent. Once a task is 
started by an agent, it will be finished completely by that 
agent at one time. Wien an agent completes a task, the reso- 
urces possessed by the agent are consumed and a certain amount 
of cost is incurred. The model is used to determine the best 
assignment of tasks to agents so as to minimise total system 
cost while satisfying the agent resource constraints. 


Mathematically, the model can be represented as follows s 


Min Z E c . . X . . 
iel deJ 

s.t. 

E X. . = 1, V d e J, (1.1) 

iel 

(p) .2 a^^ x^^ < V i e I, V k = l,...,p, (1.2) 

j s J 

Xid = (0, 1), V iel, V d e J. (1.3) 

This problem becomes similar to the Generalized Assign- 
ment Problem when the number of resources possessed by each 
agent is equal to one. 
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1,3 APPLICATIONS; . 

The MRGAP has applications in determining optimal 
assignments of tasks to agents when the number of tasks is 
larger than the number of agents and the number of resources 
possessed by each agent is more than one# For example, the 
agents may be the dock yards and the tasks to be done could 
be the damaged ships. The resource constraints coiiLd be the 
number of berths available, number of equipment hours availa- 
ble and the number of man hours available etc. The objective 
function co\ald be the minimisation of the cost of transporta- 
tion of damaged ships to the yards i.e. the ships have to be 
transported to the yards and the cost of transportation vary 
with the distance and size of the ship. A. similar application 
could be in the case of railway workshops and wagons. 

Alternatively, the agents could be the district head- 
quarters and the tasks are the schools in that area and our 
aim is to assign the schools to these headquarters, ignoring 
the geographic boundaries of the districts, such that the sum 
of weighted distances of these schools from their respe«tive 
headquarters will be minimum. The resource constraints for 
each headquarter could be the numbers of children it tan con- 
trol, number of inspector days available with it and the amount 
of funds etc. Similarly, police districting can also be fonnu- 
lated as ]yiRGAP. Here the agents are the police stations and 
the tasks are the localities and the aim is to assign the 
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localities to the police stations so as to minimize the sum 
of distances of the localities from their respective police 
stations, while taking care of resource constraints like the 
number of police, number of crimes that can be controlled and 
the amount of funds available for each station. 

Another important application of MRGAP is the optimal 
assignment of jobs to computers in a time sharing computer 
network environment. The agents are the computers and the 
tasks are the various programs that are to be executed on 
these computers. Each job may take different processing times 
on different computers and it may cost differently on different 
computers. Also, the various types of memory requirements 
may change from computer to computer, particularly when the 
computers are of different categories like super, main frame 
or mini computers etc. So the resource constraints will be 
the amount of memory available of each type and inpxit, output 
capacities etc, and the amomt of processing time available to 
each computer. The objective function could be the minimiza- 
tion of the cost of processing a given set of jobs. 

1.4 OUTLINE OF THE THESIS? 

A brief outline of the thesis is given in this section. 

In Chapter II, ve survey the existing literature in 
the field of Generalized Assignment Models, We present the 
salient features of the Ross and Soland’s algori1±im [ll] for 
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the generalized assignment problem vhich is a special case 
of 'OUT model and discuss the advantages and disadvantages of 
his approach to our model. We also give the applications of 
lagrangian relaxation from the literature. Also, a brief 
account of the concepts of lagrangian relaxation and subgra- 
dient optimisation is presented. 

In Chapter III, we develop the algorithm for MRGAP. 

The algorithm is developed on the basis of branch and bound 
approach. In this algorithm, we use the lagrangian relaxation 
concepts to a limited extent which improves the performance 
of the algorithm substantially. Another feature of the algo- 
rithm is that the lagrangian multipliers are not updated in' 
the latter part of the branch and bound tree. The search 
method adopted in the algorithm is depth first search. It has 
been found empirically that depth first search will yield a 
good feasible solution early in the branching process. 

To test the computational performance of the algorithm, 
a FORTRAN program was -written and tested on a number of randomly 
generated problems. Computational results showed that the 
algorithm can solve' reasonably large sized problems. It solves 
a 1000 0-1 variables and 6 resovirces problem in about 60 seconds^ 

These computational results are presented in detedl. A numeri- 
cal example was also given in this chapter. 

Computational results have also been presented f oif. a 
e - optimal algorithm -vdiich is useful for solving -very large 
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sized problems i It was found that with a value of 0,01 for 
e, the algorithm has been able to solve a 8000 0-1 variables 

and 6 resources problem in about 3 minutes. 



CHAPTER II 


LITERATURE REVIEW 

2.1 INTRODaCTIONs 

In this chapter, we present the previous work done in 
the area of the assignment models. We also present, in detail, 
the Ross and Soland’s algorithm for GAP and discuss its 
merits and demerits to apply their approach to our model. 

We then give a "brief review of literature on the lagrangian 
relaxation and subgradient optimisation* We also briefly 
explain the concepts of lagrangian relaxation and subgradient 
optimisation at the end. 

2.2 LITERATURE REVIEWS 

The development of algorithms for various assignment 
models has been a potential area for research for a long time 
and continues to be one. In 1956, KUHN [lO] developed an algo- 
rithm for the classical assignment model which is known as the 
Hungarian method. Subsequently, several alternate algorithms 
have been suggested for this problem. Barr, Glover and 
Klingman [2] have designed an alternating path basis algorithm 
to this model which does not consider all feasible bases for 
progressing to an optimal basis. 
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Hung and Rom [s] suggested an algorithm based on a 
scheme of relaxing the given problem into a series of simple 
network flow problems for each of which an optimal solution 
is easily obtained. 


Another model, namely Generalized Assignment Problem (GAP),: 
which is a generalization of the classical assignment problem, 
has been formulated by Ross and Soland [ll] in 1975. The 
mathematical formulation of this model is as follows s 


Min S 2 c. . X. . 
iel jeJ ^ ^ 


s,t. 

2 X . =1, V j e J (2.1) 

i£l 

(GAP) 2 a.. X.. < b,., V i e I (2*2) 

jeJ 10 10 - 1 

x^^ = (0,1), V i e I, ¥ 3 e J. (2.3) 

This problem takes on more of the appearanoe of the 
classical assignment probl^ when each constraint (2.2) is 
scaled by dividing both sides of the inequality by the right- 
hand side value to obtain, 


2 K, . X.. < 1, ¥ i e I (2*4) 

jeJ 10 10 

where, K. ^ = Si.j./'b. 

1 J 1 J X 


It is apparent that the classical assignment p3X)blem 
is a special case of GAP in which * 1 for all iel and 
e J and fa = n. 





11 


Ross and Soland have siaggested a branch and bound 
based exact algorithm for GAP, The main featin?es of their 
algorithm are that the algorithm solves a series of continu- 
ous value knapsack problems in order to obtain tight lower 
bounds. The validity of this lower bound has been piroVed 
by viewing it as a specific application of lagrangian relaxa- 
tion where the constraints (2.1) in GAP are relaxed. The 
objective function value of this relaxed version will provide 
a lower bound for the optimal solution of the original prob- 
lem GAP, This approach of lagrangian relaxation is commonly 
used for constructing tight lower bomds in integer program- 
ming problems. 

Another characteristic feature >hich has beaa exploited 

in their algorithm is the LIFO treatmoit of candidate problems 

in the branch and bound tree. This is beneficial both in tenns 

of core storage and solution time requirements. With this 

arrangement only those restrictions associated with the current 

* 

relaxation need to be stored. It has been observed by them 
empirically that LIFO approach yields a good feasible solution 
early in the branchirig process. 

Later, other algorithms have been suggested for GAP. 
Fisher, Jai Kumar and Van Va^senhove' [4] have developed a 
multiplier adjustment method. and Klastorin [9] developed a 
hetJristic method for the same problem to solve large sized 
problems. 
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A general formulation of these models have been given 
by Ross and Zoltners [12] who gave a general name weighted 
Assignment Models to this type of problems* The mathematical 
formulation of this general model is as followss 


Min r Z S c . X. -v. 
iel jeJ heH. . 


s.t* 


(WAP) 


Z S ^iih * ^ d £ 

iel heH.. 

3? . ^ < S 2 . X. < 

^ " jeJ beH^^ ^ * 


(2.4) 


V i e I, ¥ k e K, 


(2,5) 


Xijh “ ^ iel, VJeJ, ¥he (2,6) 


This model represents problems with the following charac- 
teristics: A set of tasks must be divided among a set of agents, 
and each task must be completed by only one agent at one of 
the several predetermined levels. When an agent completes a 
task at some level, the resources possessed by the agent are 
consmed and a certain amoimt of cost is incurred. The model 
is used to determine the b^st assignment of tasks to agents 
and the task completion levels. 

It is clear that the classical assignment problem, GAP 
and MRGAP are the special ^sases of WAP, When the tasks are 
Completed at . only one level, that is, if a task is assigned to 
some agent, it will be finished completely by that agent at a , 
time and when the number of resources associated with each agent 
is eqiml to one, the resulting problem will be GAP# 
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Ross and. Zoltners also discussed the relation of this 
model to the transportation models, knapsack models, facility 
location models, transignment models etc. They have provided 
applications for various models described by them and the 
information about their algorithms. 

The model considered in ttiis thesis (MRGAP) is also a 
special case of WAP where the tasks will be completed at 
only one level. The number of resources associated with 
each agent can be more than one. 

The difficulty in applying the approach of Ross and 

Roland to our model is that the knapsack problems will become 

\ 

multi-dimensional in oixr problem. Each knapsack problem will * 
have p constraints where p is the number of resources asso- 
ciated with each agent. We are not aware of any algorithm 
for the multi-dimensional knapsack problem. 

Because of this difficulty, a different approach of 
the lagrangian relaxation has been adopted in this thesis -vdiere 
instead of constraints (l.l), aS has been .done, by Ross and 
Soland, we relax constraints (1.2) in (P) , The lagrangian 
relaxation concepts and their applications in literature have 
been excellently reviewed by Fisher [3] in 1981, He gave details 
of the basic .constructions of lagrangian relaxation and its 
applications. He also" discussed the criteria to decide between 
the competing relaxations as thej^ arfe- more than biie wsCy of 
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applying lagrangian relaxation to a model* He also reviewed 
the merits and demerits of various methods of updating the 
lagrangian multipliers and finally explained how these con- 
cepts can be used in a branch and bound algorithm* 

The lagrangian relaxation has been \ised successfully 
to various models in the literature. Held and Kacp [6] have 
first applied this principle in their highly successful branch 
and bound algorithm for travelling salesman problem. Subse- 
qiaent to the successful work of Held and Karp, many applications 
of lagrangian relaxation have appeared in the literature. 

Shogan [13] applied these concepts to the resource constrained, 
capacitated minimum cost spanning tree problem, Gavish [5] 
applied lagrangian relaxation to a computer network design 
problem. 

The updating of lagrangian multipliers is the critical 
part of the lagrangian relaxation procedure. Obtaining the 
optimal multipliers, that is, the one which gives the highest 
lower bound, is highly time consiraiing owing to large number of 
linear inequalities involved. So a new method, namely subgra- 
dient optimisation is being used widely in the literature with 
tremendous success. In this method, the aim is to. calculate 
a suboptimal set of multipliers, by approximating the sub- 
gradient as the steepest ascent gradient to the optimal solur- 
tion, A detailed discussion of subgradient optimisation is 
given by Held, Wolfe, Crowder [7], Bie proof of itxe validity 



15 


of the subgradient optimisation procedure is also given in 
this paper, 

2.3 LAGRANGIAN RELAXATION CONCEPTS? 

In this section, we will briefly e3q)lain the concepts 
of lagrangian relaxation and its validity. 

We begin with a combinatorial problem formulated as 
the integer programs 

Z = min c x 
(PL) Ax < b 

X > 0 and integral 

We assimie that the constraints of (PL) have been relaxed 
to make it easy to solve the lagrangian problem 

Zq (v) = min (cx + v (Ax - b)) 

(LR^) X ^0 and integral 

where, v is the set of lagrangian multipliers. 

Now, we show that the objective function value of (iJRy.) 
will be a lower bomd on the optimal value of Z for any parti- 
cular set of lagrangian multipliers, i.e, we want to show that 
Zj^(v) <, Z for any v. 

When the constraints of the type Ax < b are relaxed, we 
reqiiire that v ^ 0. 

Now, let x*be the optimal solution to (PL), then 
Zjj(v) < (c X* + V (Ax* - b)) < Z 
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where the first inequality follows from the definition of 
Zq(v) and the second inequality from Z = c x*, v ^ 0 and 
Ax* - b ^ 0, Similarly, for Ax ^ b, we require v ^ 0 for 
Zj^(v) ^ Z to hold. 

The fact that Zj^(v) ^ Z allows (LR^) to be used in 
place of (LP) to provide lower bounds in a branch and bound 
algorithm for (PL) . 

Determining v ? 

As a tighter lower bound will help the branch and bound 
algorithm in pruning the nodes effectively, our aim is to 
determine v such that Zj^(v) will be maximum. In other words, 
our aim is to solve the problem, 

Zq = max Zq(v) 

(D) ^ 

V > 0 

There are several methods to determine the optimal v, 
but according, to Fisher [3], the subgradient method is the 
simple to program and also requires minimum computations and 
gives very good values of v. We briefly explain the subgra- 
dient method here. 

In the subgradient method, given an initial set v°, 
a sequence {v^} is generated by the rule 

(A:c^ - -b) 
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•where is an optimal solution to (IRy) and t^ is a positive 
scalar step size. 

The step size used most commonly in practice is 

t = Xr - Zd(v|^)) 
ll Ax^- bll^ 

where is & scalar satisfying 0 < < 2 and Z* is an 

upper bound on Zj^, 

Unless we obtain a v^ for which ^(v^) equals the cost 
of a known feasible sol'ution, there is no way of proving opti- 
mality in the subgradient method. To resolve this difficulty, 
the method is usually terminated upon reaching an arbitrary 
iteration limit. 



CHAPTER III 


THE BRANCH AND BOUND AIDORITHM 

3.1 INTRODUCTIONS 

In this chapter, we develop the branch and bound algo- 
rithm for MRGAP based on the lagrangian relaxation concepts. 
Some of the important features of the algorithm are that 
(l) the lagrangian multipliers are calculated only at the 
root node, (2) the lower bounds ^e not calculated at all the 
nodes of the branch and bound trae. 

In Section (3.2) we present the flow chart of the 
proposed algorithm and a step-by-step explanation of the flow 
chart. In Section (3.3), a numerical example is given. In 
Section (3.4) computational performance of the algorithm is 
presented in detail. The results of e-optimal procedure are 
also presented. Finally, concluding remarks are given in 
Section (3.5), 

3.2 THE BRANCH AND BOUND ALGORITHMS 

The branch and bound algorithm developed in this thesis 
has a special feature that the bounds are calciilated using 
lagrangian relaxation concept $ and some additional characteriS' 
tics have been exploited to reduce the computations at each 
node , 
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The basic idea behind any branch and boimd algorithm 
is to investigate'- all the feasible solutions by partitioning 
the solution set into subsets, by a process known as branching, 
A lower boimd is calculated for the ob;)ective function value 
of the solutions in each subset. If the lower bound is greater 
than or equal to the ob;i active function value of the incumbent 
solution, that is, all the solutions in the subset are worse 
than the incumbent solution then further branching from this 
subset is stopped. This is called pruning. The partitioning 
continues until the objective function value of the incimibeirb 
soliation is lower than any bound for any subset. The algorithm 
terminates when all the subsets are primed. 

In the discussion to follow about the branch and bound 
algorithm, Z denotes the objective function value of the incum- 
bent, the best known feasible solution of (P), Z is initially 
set to the objective function value of the solution obtained 
by assigning each task to the agent with maximum cost ignoring 
the resource constraints. 

In the branch and bound tree, the root node represents 
the problem (P) , We will partition the solution set of this 
problem over a separation variable, let us say x^^ into two 
subsets, one subset containing all the solutions with x „ 

vL * 

equal to one and the other containing solutions with x^^ equal 
to zero. These resulting problems of petrtitioh are similar to 

■ , ; ”>r 

(P) but additional restrictions are placed oh the solution 
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matrix x. These problems are of reduced size and are called 
candidate problems. Now, our aim is to calculate the lower 

ju 

bounds for each of these two candidate problems. Let (P“) 
represents a candidate problem. 

The basic concept in the bounding procedure of the 
algorithm is best illustrated with reference to the initial 
relaxation of problem (P). There are two natural relaxations 
of (P), one is by relaxing the constraints (l.l) and the 
other by relaxing constraints (1.2) in (P) . The mathematical 
formiJlations of these relaxations are given belows 


"When constraints (l.l) are relaxed, the model becomes. 


Z^(u) = Min I E x-^ + £ u^ (l - E x^. ^) 


4-1 ZJ ^.5 4 -^4 -1 

iel ;ieJ jeJ ^ 


iel 


Xj 


or, equivalently, 

Z^(u) = Min 

s.t. 


(LR^) 


I Z - uj X . + Z u 

iel jeJ J “ jeJ ^ 


^id ^Id ^ ^ i e I, ^ k = 1, 


OeJ 

Xy = 0 or 1, V d e J, i e I. 


Thus (LR^) separates into a series of m\ilti-dimensional 
knapsack problems, one for each iel. This is a hard problem 
if an exact solution is desired. However, a lower bound can 
be obtained by solving m linear programming problems of 
size p. X n. 
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The second relaxation (LR2) is obtained by relaxing 
constraints (l*2) with v ^ 0 

Z^Cv). = Min S E c.. X. . + Z E v?( Z af .x..-b?) 

iel deJ ^ ^ iel keK ^ jeJ ^ 

or, equivalently, 

^(v) - Min S E (c, .+ E vfa^Jx.,- E E v? 

jeJ iel keK ^ ^ 


iel keK 


s»t. 

^^2^ E X. - = 1, ¥ d e J 

iel 

Xid = 0 or 1, ¥ d e J, ¥ i e I 

Now, one may discuss the merits and demerits of these relaxa- 
tions* To evaluate a relaxation, we have to consider two pro- 
perties J The * sharpness of the bounds and the amount of compu- 
tations required to obtain these bo\jnds. Usually, selecting a 
relaxation involves a tradeoff between these two properties, 

It is generally difficult to determine whether a relaxation 
with sharper bounds but greater computational time requirements 
will result in a branch and bound algorithm with better overall 
performance. In our problem, the relaxation (U^) needs solving 
m multidimensional knapsack problems at each iteration or 
alternatively, m linear programming problems. A multidimensional 
knapsack problem [l4] is one with the same basic construction 
as the well known knapsack model but have more than one cons- 
traint . This problem is a well known NP-complete problem and 
we are hot aware of any algorithm for solving it. Also, it may 
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not "be computationally feasible to solve m linear programming 
problems at each iteration. 


In contrast, solving (LR2) requires solving problem of 


the type 


Min (c.. + E v.^ a^. ) x. . 



To solve this, we have to find an agent among I, which 

p k If . 

is having minimum (c. . + E vf af.) value for each e J and 

k=l ^ 

the corresponding decision variable is made equal to one. This 
requires less computations, of the order of mnp, when compared 
to (LRj_) . Because of this great difference in computational 
time requirements we adopt relaxation (LR2) in our branch and 
bound algorithm, though it produces feeble lower bounds rela- 
tive to 

Now, we will see how the lagrangian relaxation principles 
increase the efficiency of the branch and bound algorithm. 

Let 2^ ^ 2 ^^^ denote the optimal objective values 

for problems (P) and (LR2)» respectively. Clearly, Z2(v) - ^opt» 
that is, Z2(v) is a lower bound on for any non-negative v, 

eus' shown in Chapter II, A traditional branch and bound approach 
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might xise as a lower bound for that is, compute a 

lower bound on Z^^^, that is, compute a lower bound on Z^^^, 
by simply ignoring the resource constraints and solving the 
resulting problem. However, a branch and bound algorithm 
based on L-relaxation seeks to increase the likelihood of 
fathoming a node by searching for the greatest possible lower 
bound, that is, the objective will be to solve 


(PM) 


= maximise Z 2 (v) 
subject to 


vj > 0, V i e I, V k = 1, ..., p. 

There are other benefits of this approach in addition to 
the fathoming of nodes on the basis of lower bounds. In partici*- 
lar, instead of simply solving (P^), the attempt to solve (PM) 
Involves the solution of a series of problems of the form (LR 2 ) ' 
with a different value of v each time. This increauses the 
likelihood of obtaining a feasible solution which may be better 
than the incumbent value, thus increasing the likelihood of 
fathoming the nodes in the subsequent part of the branch and 
bound tree. 

Fig. ,1 contains the flow chart of the branch and bound 

algorithm. . A more detailed discussion of 4ach of tte steps 

follows s Let denote an m x p matrix whose is the value 
/ k k\ 

( Z a. . X. j ~ b.) evaluated at v. A scft>-problem In the brsoich 

jej Id la 1 

and bound algorithm consists of three sets of variables J 
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variables whose values are fixed to zero, variables i^diose values 
are fixed to one and the variables which are unassigned. For 
siniplicity of presentation, let us assume that two matrices 
y° and y^ are associated with each sub-problem defined as 
be lows 


v° - 

1, 

if Xj^j is 

fixed to zero 


0, 

otherwise, 

V i e I, ¥ 0 

-Id = 

1, 

0, 

if x^j is 
otherwise , 

fixed to one 

¥ i E I, ¥ 0 


These are used to represent the additional restrictions on the 
candidate problem at any instant of the branch and bound tree. 

Step li In the branch and bound tree the first node, i.e. the 
root node, represents the initial relaxation of the problem (P) . 
The lagrangian multiplier set is calculated at this node using 
subgradient optimisation and will be used throughout the tree. 

To start with, all the lagrangian multipliers are made 
equal to zero. In lagrangian relaxation, we try to increase 
the cost of allocating the task to those agents which are over- 
loaded, they making it possible to transfer some of the tasks 
to other agents which are under utilised. ¥e do this by adding 
a factor to the actual cost of allocating a task. This factor is 
as seen in the formulation (LR^), where v^ is the 

k*l ^ 

lagrangian multiplier corresponding to the constraint on k-th 
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resoijrce of i-th agent. It is natural to start with a value 
of zero to all v^*s and then increasing v’s of those agents 
which are overloaded. 

Now, our aim is to obtain a good suboptimal, if not 
optimal, values for the lagrangian multipliers using sub- 
gradient optimisation. For this, we make an attempt to solve 
(PM). But solving (PM) optimally is impractical owing to the 
large number of linear inequalities represented in the problem 
(LR2), hence an attempt to solve (PM) will be made throijgh an 
iterative procedure known as subgradient optimization [7], sub- 
gradient optimisation is so named because is a subgradient 
of the objective function for (PM) at the point v, and the 
procedijre optimally loses this subgradient as if it were a 
gradient pointing in the direction of steepest ascent. Itera- 
tion r of the subgradient optimisation consists of the following 
steps (1) movement to the new v’^ obtained recursively from 
v^-l ^2) solution of the L-relaxation using new v ( 3 ), deciding 
whether to perform another iteration or not. The hope is that 
the sequence (v^) will yield a good suboptim^, solution to (PM). 
Each step of the suhgradient procedure is estplained below in’ 
detail. 

Given the lagrangian multiplier set, the relaxed problem 
(LR2) is solved by assigning eaoh taskto that agent with minimum 

P V k\ ' 

(c, . + I Sj. V.) . value and whose, capacities are not less than 

k=l ^ 

the consumption of the resources by this task. 
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Af-ter (LR2) is solved with a particular set of v, Vre 
try to revise the values of lagrangian multipliers to obtain 
a bigger objective function value for (LR2) which is a lower 
bomd on the objective function value of (P). Updating lagran- 
gian multipliers is the most complicated procedure in the 
branch and bound algorithm. The procedure is based on a 
proof [7] that if the non-negative scalar t^ is sufficiently 
small, then 

v^"^^ = max [0, (v^ + t^ s’'^)] 

is closer than v^ to the optimal solution of (PM), although 
it does not guarantee an improvement in the lower bound. The 
step size t^ most commonly, used in practice is, 

_ (Z* - 2g(v’^)) 

I |Ax^ - -bl 1^ 

where is a scalar satisfying 0 < < 2 and Z* is ah 

upper bound on the optimal solution value of the problem (P), 

The value of can be obtained empirically with some esiperi- 
mentation and it is halved whenever Z2 has failed to improve 
in some fixed number of iterations. 

Now, a decision should be taken as to when ths subg^'adient 
optimisation procedure is stopped. There is no way of proving 
the optimality of v»s iii the subgradient method unless we obtain 
a V such that Z2(v) is equal to a known feasible solution value. 
To resolve this difficulty thS method is usually t^pninated 
iq>on satisfying one of the fallowing conditions. 



SOLVE INJTIAL 

LAGRANGIAN RELAXATION PROBLEM 
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(a) The scalar step size 

(b) The scalar < £2 

where and £2 are very small quantities or, 

(c) the number of iterations r > D, where D is an arbitrary 
iteration limit. 

Step 2.3.7 ; The stack contains "the list of unfathomed prob- 
lems. A subproblem in the stack is, rep resented by a vector of 
information (y^,y°,v, 2^,3"''^) , where v is not necessarily the 
optimal solution to problem (PM) but is such that Z*^ is. the 
greatest known lower bound on the optimal solution value of 
the subproblem (P ). 

Initially, in Step (l),. we solve, the problem (LR 2 ), at 
the root node of the branch and boxmd tree, which results from 
ignoring the resource constraints in problem (P)» Then the 
partitioned subproblems from this root node are stored in the 
stack. The storing of subproblems is done in such a way. that 
it needs imi niniiiTn effort "vdien a variable is selected to branch 
out from a particular node. 

In selecting a new subproblem from the stack, the 
subproblem is chosen according to the depth first r\Jle, that 
is, the chosen subproblem is the one which is the latest 
addition to the stack, or equivalently, the one on the top 
of the stack. This approach results in a significant reduction 
in the computer storage, since only the details of the latest 
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sufc^roblem have to be stored at any instant of the algorithm. 
However, if the stack is empty (Step 7), that means all sub- 
problems have been pruned and the incumbent solution matrix 
is the optimal solution of (P) , 

There are other strategies existing for selecting a 
subproblem from the stack. The most appealing alternative is 
the best-bo'und rule, which selects the new subproblem as the 
one with the lowest lower bound among the unfathomed subproblems. 
Such an approach would result in lower percentage errors* when 
the algorithm is forced to terminate prematurely due to excessive 
execution times, because we are considering the most promising 
subproblems first. Ihsp5.te of this, the depth first rule was 
chosen because the best-bound rule requires lot of computations 
to i^date the list of subproblems and also because it is observed 
empirically that the deplh first rule will yield a good feasible 
solution early in the branching process. Of course, there is 
no guarantee that the depth first rule is the best choice. 

Step 4; (Computing the lower boxmd) An important feature 
of our branch and bound algorithm is that the lagrangian 
muiltipliers are not updated at each execution of this step. 

This approach is motivated by the intuiitive consideration that 
the computations involved in determining optimal m\£Ltipliers 
are too time cpnsuuning and may not ^uistify the improvement 
obtained in the lower bound value. ,We thus. prefer to sacrifice 
the sharpness of the bound to achieve the speed of executing 
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this step. This approach is also found to be highly conducive 
to reoptimisation from one subproblem to the other. Since 
Vp once confuted initially, are fixed throughout, the ( 1 ^ 2 ^ 
reduces to 


Zo « min 2 Z C . . x . . 
2 OeJ iel 

sub;3®c'fc "to {1*1) and (1,3) 


2 2 
iel keK 



where. 


* keK 


The relaxed subproblem to be solved is 


Min 2 2 c . y}. + 2 2 C. . x. . - 2 ^ ^ 

iel jeJ iel jeJ iel keK 

subject to (l.l) and (1,3) 


where, 




5^-2 a^^ y}., V iel, V keK 

'i xj ^10' 


Clearly, the lower bound for the subproblem is 

LB = 2 2 c . . yL - 2 2 • b^ 

iel jeJ iel keK ^ 




jeJ, 2 y^ =0 
iel 


r 

min 

iel and 
j^and (i, j);^ S 

.k 


where S is the set of variables whose a;^ > b^ for any k. 

Now, from the above expression it is clear that when 
a variable is fixed to one or zero, the computations needed to 
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update the terms in the .expression are very less. For example, 
when a variable x^^ is fixed to one, the second term in the 
expression is updated since b^ changes. Then it is checked 
whether there is any task j, which is assigned to agent a pre- 
viously, vrfiose a^ . value is greater than the new b^ value for 

U J Wr 

any k. If there are suchtasks they will be reallocated and the 
lower bound will be x:ipdated accordingly. This needs computa- 
tions of the order of qm, where q is the number of tasks to 
be reallocated. Otherwise the lower bound will not change and 
so it need not be recalculated, which is the case when a fea- 
sible solution is found. Similarly, when a variable x^^ is 
fixed to zero, the second term is updated and the taskwill be 
reallocated and the third term in the expression is updated 
accordingly. Also, we check whether any tasks have to be re- 
allocated to this agent because its b^ value has increased 
when we fixed the variable x^^ to zero. Thus, the lower bound 
need not be calciiLated afresh at each node> but it needs to be 


updated only. 


Step 5S After solving (P^) , we will try to prune the sub- 
problem that is, to remove the subproblem from further consi- 
deration. This pruning can be divided into four steps. In 
the first step, we make a check whether ^ , the lower bound. 


is greater than Z, the incumbent solution value at that instant 
of the branch, and bound tree. If it is so, the. subproblem is 


pruned ^ it is removed from further consideration. However, 


I tv . 


■'TTro 
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if "tliis condition is not satisfied, we proceed to the second, 
step* In this step, we check whether the solution x(v) is 
feasible to (P), that is, whether the solution x(v) satisfies 
the resource constraints of (P) also. The feasibility check' 
can be made easily by checking for a' strictly positive com- 
ponent in S^. If all S'^ < 0 does not hold, the pruning 
process ends without pruning the subproblem. If all s"''’ < 0, 
we proceed to the third step, where we check whether the 
solution of (P^) is better than the incumbent solijtion. If 
it is so, the incumbent solution is updated, 'Whatever is the 
result in the third step, we proceed to the fourth step. Here, 
we will make another atten^t to fathom the subproblem. When 
x(v) is feasible to problem (P) also, the following condition 
will be satisfied. 

« cx(v) + vS^ <_ <. cx(v) 

where the last inequality follows from the verifications in 
the second step of the pruning process, of the feasibility 
of x(v), Also all ;< 0, because x(v) is feasible to (P). 
Hence, if = 0, then 

^ ■= ^opt = 

so that V and x(v) are optimal solutions to (PM) and the sub- 
problem (P^)» respectively* Note that, if subproblem (P^) is 
fathomed in the fourth step of pruning process, then x(v) will 
al'^ys have become the new incumbent, because 
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cx(v) 2?^ < 2 

whe3?e ■fche first: equali'fcy follows because 2^ = ^opt: ~ cx(v') and 
"the second inequality follows "because the subprohlem was not 
f Ethoffi^d in “the fips*t step of the pruning process • 

S-fcep 6; After the pruning attempt failed, (P^) must be parti- 
tioned into two subproblems. For this, we have to choose a 
separation variable x^^ to partition the solution set of (P^) 
into two subsets. There are several strategies to select a,p 
and a few are explained below. 

ha. the strategy which we adopted in our thesis, a, p 
are selected as follows, a is the one among the set of agents 
for vtoich 

T « max [ 2 max [O, 2 a^. x. . - b^]] 

® 1 k«l deJ ^ 

In other words, we are selecting the agent with maximum 
total exceeded capacity. 

Then p is selected from among the tasks that are assigned 
to a a nd which has the maximum penalty p^. The penalty p^ is 
defined as the excess cost that will be incurred when the task 
J Is shifted from the present agent to the agent with the next 
highest cost. 

The variable chosen by this rule represents a task a 
that is best kept with agent i considering both the penalty 
for switching the tasks and the resources available to the 
agent. 
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There are other strategies for selecting the separation 
variable. One such strategy is to select that agent whose 
resource has exceeded by a maximum amount than any other 
resource of any agent. Then selecting the task is same as the 
above strategy, that is, the highest penalty task. This stra- 
tegy also worked well though not as good as the first strategy. 
This is because an agent, all resources of which has exceeded 
the os^acity, is ignored, 

3,3 NUMERICAL EXAMPLES 


In this section we give a small numerical example and 
the branch and bound tree for the same. The objective function 
values (c^j)aM the coefficients (a^^) of the multiple assign- 
ment constraints for a problem with three agents, five tasks 
and two resources are shovai in Tables 1 and 2 respectively. 

In this example, the right hand side values (b^) of the 
multiple assignment constraints are 38 for all i e I and k e K, 


Initially (LR 2 ) is solved to obtain the solution x^^. “ 

Xj^g » 1 , « 1, X 24 = 1, X 25 = V and all other x^^ = 0. This 

solution is not feasible to (P) . The lower bound (LB) provided 
by this solution is 100 . 

» ^ 

The variable Xg^ is selected as the separation variable 
and the candidate problem with X 2 ^ = 1 is solved next because 
there are no other tasks which are assigned to agent 2 and which 
can be fixed to the same agent. The solution to this candidate 
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problem is = 1, ^^3 = 1» ^ 2 ^ = li ^^5 = 1 and 

all other = 0, vdiich Is again not feasible to (P)* The 
LB is 102. 

Now, variable x^^ is selected for braiiching and the 
candidate problem is solved to obtain the solution x^^ - 

^12 * ^13 ^24 “ ^55 ” ^ other x^^ = 0 

with LB « 106. 

Now, the variable x^^ and then x ^p are , selected as sepa- 
ration variables ani are fixed to one and the next candidate 
problem with x2/j, * If = If x^^ = 1, x ^2 “I* solved at 
node (5) to obtain the solution x^^= 1, x^2 “ ^33 ^*^4 “ ^ 

Xj^ w 1 and all other x^^ = 0, This solution is feasible to (P) 
Table H Objective function coefficients 







Tasks 




i*”* 



3 


4 

5 


1 


22 

16 


37 


45 

34 

Agents 

2 


32 

19 


47 


22 

15 


3 


11 

42 


34 


43 

15 

Table 

! 2: 

Muiltiple assignment coefficients. 





k«l 





k=2 





Tasks 





Tasks 


1 

2 

3 

4 

5 

1 

2 

3 

4 5 

1 

13 

6 

10 

6 

7 

8 

18 

15 

18 22 

Agents 2 

19 

13 

6 

7 

12 

12 

21 

15 

23 21 

3 

17 

24 

21 

21 

10 

25 

15 

16 

19 17 
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with an ohjective function value of 109 and the inoumbent is 
xapdated accordingly* Now, we fix all the remaining tasks which 
are not fixed in the previoixs candidate problem and back track. 
While back tracking the candidate problem with.x^^ = 0 and all 
the remaining tasks fixed is considered next. This gives a 
lower bound of 115 and is prmed. The process is repeated 
until all the nodes are pruned and finally the optimum solu- 
tion is ** "^2 ” ^33 * ^2 A- ~ ^35 ” ^ and all 

other x^j « 0, The solution value is 109. In this case, the 

first feasible solution is also the optimum solution. 

The complete branch and bound tree for this example is 
given in Fig. 2. The numbers in the circles indicate the node 
numbers thB designates the lower bound on the node and Z is the 
incumbent value. 

3.4 COMPUTATIONAL RESULTSi 

In this section we report the computational performance 
of the branch and bound algorithm for MRGAP in solving problems 
of different sizes. We also report the performance of the 
e - optimal algorithm with e = 0.01. 

The branch and bound, algorithm was coded in FORTRAN IV 
and tested on DEC 10 computer system under multi programming 
envircoiment . The program was executed on randomly generated 
problems with yaryihg number of agents (m), tasks (n) and 
resoin:*ces (p). .We considered six combination of tasks 
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sgeniis for eacli coiribina.1:ioii "three problems were solved 
with the number of resources equal to 2, 4 6* The prob- 

lems were generated using the same scheme as used by Ross and. 
Soland to generate generalized assignment problems. The 
values of and a^^^ were generated from a uniform probabildty 
distribution with intervals as [10,50] and [5, 25] respectively * 
To establish binding multiple assignment constraints, each 
b^ was set equal to 0.6 x ~ x 15 + 0,4 x R, -vdiere 


R m max [ Z a^ . x. .] 
iel deJ 


and are the solution to the initial relaxation of (P) , 
Ross and Soland reported that such a scheme attempts to set 


tight constraints and avoids the possibility of generating 


trivial problems which might be solved in the initial relaxa— 

tion. m general, based on the range of .values of the 

one would expect random problems to be infeasible if each 

b? < - X 15 and trivial if each b^ 

i m 1 "• 


For each problem solved we noted the number of nodes 
enumerated in the branch and bound tree and the conputational 
time • The times given in the tables wei*e measured by real 
time clock, accurate i;pto one millisecond, and does not inclTjude 
irput and output times. We also noted the number of nodes 
examdlned ajxl the con^jutational time taken and the objective 
function value for the first feasible solution. All these 
details are given in Table 5. Wie performance of the e-optimal 
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algorl'thin (the "branch and bound algoid-bhin in which a node is 
pruned whenever LB x (l + e) ^ incumberrf: value) is given in 
Table 4# The following conclusions can be drawn from bhe 
tables. 

1« The algorithm can solve problems U 5 )to 1000 variables 
and 6 resources within reasonable computational effort* The 
computational times are found to be more sensitive to the 
number of agents than the number of tasks. Further, the 
con^jutational times do not appear to grow exponentially with 
increase in the number of resources. In fact, increasing the 
number of resources reduces the number of feasible solutions 
and, thereby, reducing the amount of in^licit search. 

2* The algorithm examines a fairly large number of nodes in 
the branch and bound tree before terminating. This is because 
the lower bound provided by the relaxation (LR 2 ) is somewhat 
loose. However, the bound is obtained with little computa- 
tional effort which makes up the loosen^^s of the bound* 

3. The algorithm obtains a very goO^ feasible solution, 
which is often optimum, early in the implicit search. The 
optimum solution is also invariably obtained in the early 
stages and rest of the time is spent in proving the optimality 
of the solution. Hence, premature termination of the algorithm 
can be recommended for large sl^ed problems. 
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4* The performance of "the e *- op'timal algorithm is also 
very encouraging. The algorithm is able to solve problems 
i^to 8000 variables and 6 resources ■within 1 percent of the 
optimality in less than 3 minutes* 

3.5 CONCLUBING REMARKS? 

In Chapter I, ■we observed that MRGAP has useful practical 
applications. We gave the statement of the problem and the 
applications , 

In Chapter II, we revie'wed the literature available 
for the assignment models and the lagrangian relaxation con- 
cepts we also gave a brief accoimt of concepts of lagrangian 
relaxation and subgradient optimisation. 

We developed a branch and bound algorithm in Chapter III 
with special features like calculation of lagrangian multipliers 
only at the root node and calculation of lower bound only at 
intermediate nodes etc. 

Finally, we have done extensive computational study on 
the algorithm. We observed that the algorithm is performing 
quite well. It was able to solve 1000 0-1 variables problem 
with 6 resources in about 60 seconds of CPU time. We also 
observed that the e - optimal algorithm with e = 0,01 was very 
effective in solving large sized problems, and it was able to 
solve a 8CX)0 0-1 variables problem with 6 resources in about 
3 mdmites. 
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Taible 3* Performance of the exact branch and bound algorithm 
\ computational times are in seconds on DEC 10)* 


Prob- First Feasible Solution Optimum Solution 


lem 

No. 

m 

n 

P 


Nodes 

Exami- 

ned 

Compu- 

tational 

Time 

Node's 

Examined 

Compijt at ionsil 
Time 

1 

5 

100 

2 

1.000 

0 

0.39 

495 

0.83 

2 



4 

1.001 

44 

0.77 

2179 

5.30 

3 



6 

1.012 

80 

1.09 

8042 

13.86 

4 

10 

50 

2 

1.000 

25 

0.44 

1258 

2.61 

5 



4 

1,022 

37 

0.75 

176 

1.18 

6 



6 

1.000 

37 

1.06 

2005 

• 3.86 

7 

10 

100 

2 

1.001 

11 

0.84 

113 

1.15 ■ 

8 



4 

1,000 

30 

1.46 

160 

1.97 

9 



6 

1.000 

99 . 

2.09 

135 

2.68 

10 

15 

30 

2 

1.009 

23 

0.57 

384 

1.33 

11 



4 

1.040 

22 

0.99 

1210 

3.84 

12 



6 

1.000 

9 

1.37 

3001 

8.39 

13 

15 

75 

2 

1,000 

26 

1.37 

11270 

20.04 

14 



4 

1.000 

12 

2.31 

75 

2.62 

15 



6 

1.001 

23 

3.32 

165 

4.08 

16 

20 

50 

2 

1.003 

16 

1.59 

4785 

8.21 

17 



4 

1.000 

36 

2.74 

500 

3.57 

X8 



6 

1.000 

18 

3.86 

13925 

54.72 


liable 4: 


Performance oi 
with e =s 0,01. 
(Computational 


Problem 

No. 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


10 100 

10 200 

20 100 

20 200 

40 200 


15 


2 

4 

6 

2 

4 

6 

2 

4 

6 

2 

4 

6 

2 

4 

6 
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the e - optimal algorithm 
times are in seconds on DEC 10). 


Nodes Computational 


Examined 

Time 

22 

1.14 

60 

1.90 

104 

2.65 

200 

1.98 

124 

4.42 

174 

6.08 

190 

3.47 

132 

5.89 

880 

8.30 

112 

5.25 

194 

12.36 

200 

17.45 

76 

8.36 

52 

45.99 

88 

175.67 
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THIS Pi^ 3G;RA'< GKNieRA^TES AND SOLVES MiiLTI-RES JijKZi GE 'i ERA u ISKJ 
A.SSt.3.N^iE^r PRDBLS*'1S. 

THE i.NpJr paramsiers are Explained sslo^ 

^ vu'Iber of agents 

N :;-'• NU'^BER OF r.ASiCS 

p :!-:- NUMBER OF resources 

isesO'J:— seed 

lAlpjRTANr variables IN THE PROGRAM ARE EXPLAINEO 3SL3aI 
1,233 incumbent solution value 

IX ss— incumbent solution 

A ij— multiple assignment constraint coefficient s 

:■ oejEcriVE function coefficients 

3 * s — cApAcrrv matrix 

P8 CAPACMV matrix QF THE SUBPROBLBM AT ANT INSTANT 

XCIf ji) * 1' IF task J is assigned TO AGENT I 
» 3 OTHERNIse 

re ss-:- sum of Aci,j,K)»x(if j) for all i and it in the 
initial calculations at the root node 
A(I, 0rKI4X(I, J)-p8(I,K) there after 
V n-.- LAGRAMGIAM MULTIPLIER MATRIX' , 

C1HE;*J ij-,- sum of ACI ,J,K)'l'V(IrK) FOR ALL' I AND J 
SE J{— PEMALTlf ARRAI 

KX It-:- heap to order the tasks in thE decreasing order of 
penalty 

FLCLfJO » TRUE IF VARIABLE X(I,J) IS FIXED TO ZERO 

= false otherwise 

PLJiS^J!) s. TRUE IF TASK J IS FIXED TO SOME AGENT 
». FALSE otherwise 
IS II-:-; STACK OF agent Indices 
JS tt-K STACK OF TASK INDICES 

IVAL< M-;-; POINTER TO INDICATE THE SUBPROBLEM STATUS 
C0NSi(J3M-K agent to which TASK J IS ASSIGNED 

^if;^1^$$.^*t**’l^****** **$***■**********************¥***** **¥****** 

COMMON »li,N,P,A(25,209,6),C{25,200),BC25,6),lX(2DD),I23D 

INTIGER' P • 

real* large.', I20D# max 

RE:AD'(21,M')M,N,P,ISEED 

CALL' G.SNRAT(M,NrP,A,C,B,ISEED) 

CALL! MR3AP 
e:izo3«izoO' 

NRlTS'C23rl003IlZO0 



90469 


4PK.S1 

90479 


^p =:23 

90489 


IMsM 

90499 

no 

lFCr.4.iir,20J^P=MpKtIN-l 

90599 


/jRir2:t2 3, 101 HI,I = '1PK,MP) 

90519 

101 

FORMi^TC/U, ' J3BS'r3X,20CI4,lX)) 

90529 


^«Rirs;(23no3) 

90539 

103 

FaR4 4T(/,8XnOClXf ' — %1X)) 

90549 


j<Rirs:(23, 1 04) 

90559 

104 

F3R'<Air(8X/20('i '>3xni ')) 

90569 


/«Rirs:t23n05)ClXtl) ,I=«PK,MP) 

90579 

105 

FOR'^Air(ixn4c;Ef4rs',iXf20('m3ni 'n 

90589 


/4RI rs;(23nob) 

90599 

106 

F3R4AT(8X#20C1X, IX)/) 

90699 

. 

lNaI.'I-20 

90619 


'IPKSMPfl 

90629 


HPa'<,p+20 

•0639 


lF(I.M,3.r.0)G3 TO UO 

•0649 


ST3P?E.?1 0 

90659 


SUBROUrifiS. 'IRGAP 

90669 



90679 

♦ ♦ 

THIS SU9R3UriNR IftKES THE PROBLEM DATA AN3. Gl^ES I'KS 

90689 

♦ ♦ 

opri.'itJ^ soLono'J ro a m-r-c-a problem 

90699 

f « 


90709 


ODM^OM »<,M,P,A(25»20p»6Ji-C(25,200),BC25,6)nX(25))nzOD 

90719 


OmENSIPN PBC25,6),V(25,6) ,rEt25,6) 

90729 


OTMfNSIDN KX(200),SE(200),lS(1000),JS(lOOO},IVAU'(l3 0a) 

90739 


0-IMSIDN C«EW(25f200) 

90749 


IHTESER X(25,200) f P, rOP,FlRST#C3NE(2aO) 

•0759 


UOGim' FCi*(25,200),FLRG(20o),FEASE,irRA:K 

90769 


REAl.' KL6,t.erIZ00vICDST, LARGE, MAX 

90779 

, f 



■0789 


EPS! so, 00001 



Sl»S2sO, 00001 

polo 


LIHItRsMi 

90819 


irCM).Cie..e)LlMlTRs21'M 

90829 


IS)bO ' ' ' 

90839 


00 7 IS1,M 

90849 


■ 00 7 KJ5l,P, . 

90859 

7 

P8(I,K.5s8CI,K3 

•0869 


LARS'SfslOOOO. 

90879 


IZOO’SO,’ 

90889 


PMI'la.lOOOO, 

90899 


00 9 I=l,M 



iC>9t>» 


33 9 ICsi^P 

00919 

9 

V(I#^)=0.0 

90929 

♦ ♦ 

****** ENlTIAlilSe IZQO ****** 

90939 


DO 22 J=1,M 

90949 


MAX*3. 

90959 


DO 21 1=1,^ 

90960 


>C(IfJ)s3 

90979 


IFCKI, J).l4E,MA)C)G3 TO 21 

90989 



90999 



91009 

21 

:DMr£>iu£ 

91010 


rZ33'=.U3DfCtI'^D'^, J) 

91029 

91039 

22 

co^riNus 

i 

91049 


DO n jslrfi ' i 

01059 


FLAS-CJls.fAliSE. 

9 IO &9 


:o«£i(j)=o, 1 

•1079 


DO 31 1 

91089 

31 

rod, J)S.FALSE, 1 

•1099 


rop=a j 

i 

•1199 


rlXir4LeO, 

91119 


irDiJ-so 

91129 


iiorAti«o 

•1130 


CAIjiU' RTIMStfll) 

OlHO 


P£*A?<OAsO,2dC3U^rsUPREVLBsO.?irERsO 

01150 


- 

•not 


***f CAiiCOLATlOS OF LAGRANGE MULTIPLIERS AX THE ROOT NDOE **** 

01170 

11 

iter»i:ter+ufease»,true. 

01180 


ICOSTftO,* 

01180 


L6s0,* 

01200 


Ei:<S^EN«0,' ■ 

01210 


DO 4o 

•1220' 


DO 4 D K«4,P I 

%m§‘ 


mi.fK3»0.' [ 

#i«o- 

40 

ELeMilM»£L€MES+va,K)»PBCIrK3 | 

01250 


llR'ACK;«.,f Al,SE.' 1 

912&9 


■ DO 30. ■ . j 

01279 


PM lN*iiARG-E( 

•1280 


DO 20 1*1# « 

91290 


TEMPsO,' 

01399 


:ne^ci:,j:)soarge 

01319 


DO 10 x:*i,p 



IFCAXI, J,X,),GT,P8CI,K)) GO TO 20 | 

r 

01339 

10 

‘rE'lP=.V£i;X.)»ACl,J,K)4-TEMP ' | 


9135» 


rEMPs.rsv(P4-c(i, j) 

9i3&9 


IFC rsMP.GE.PMlN) G3 JO 20 

91379 


pMi'is.rs'iP 

91380 


i'N0s:x=i 

91380 

20 

:3MriNUE 

91480 


IFtP'^n.GS.UARGe.) ITrACKs.TRUE, 

01419 

79 

:3'4fi.( J)=I^iDEX 

91420 


rc3srsi:osr+c(iMDex, j) 

91430 


L«sL»‘3 + pv5i\J 

91449 


00 2 5 iCsl,P 

01450 

26 

rECi.'10EX,<3 = rE(lNOEX,K)+A( index,, 

01450 

30 

:aNri;Mus 

91470 


iFCir«A:K)33 ro 16 

01480 


XLBS'tiO 

91490 


[,d»X,D8-ELSMEN 

91500 


IFCPREVD'B.GE.(J,0IJI»LB))G0 TO 45 

91510 


ir03XT*l 

01520 


PRE/LBsO© 

01530 

45 

504=0. 

91540 

' 

pMAXs^lOOOO. 

01550 


00 53 1*1,4 

91550 


00 53 X*1,P 

91570 


rSCI.,K)arE.(I,K)-Pfi(I,K) 

•1589 


xx«rai,K) 

91590 


StlH^BUX+XXfXX 

91609 


rf(XX,G,r,0)FEASE».FAtiSE, 

01619 

50 

aDuriNUE; 

91620 


IFCFEASEO GO TO 49 

91630 


ZrCWER.GT.OlMirRJGD TO 49 

01640 


rl,«F4ft40ft* C IZOD-PREVliB) /SUM 

91659 


00 53 

91660 


00 60 K.».1,P 

01670 


n i:> KD ( 1, K.) + T £ 1 * r E C I , K ) 


60 

iF(/(i;,K;).t,>r,0) vu^ioso. 

miH: 


rraw»ic30MT+i - 

91700 


IFCICDUNT.t,E,liI4)G0 TO 70 

01719 


pl,A40ftiBPI,fAMOA/2 

01720 

70 

lF(Pl4640A.6I,epSUGO TO 49 

91730 


iFcry,i,(r,sps 2 ) so to 49 

91740 


GO ro 11 

. 91750 

49 

DO 35 J;sl,S 

91769 


nosj]=:3[iE;c I) 

P 91770 


XCI-DNI, J>al 

91780 


SEGO'IO'SOA.RGE 


9l79» 


DD 85 lai,4 

91999 


IFC i:0»iJ.,23.I}G3 ro 86 

9181® 


rE'<P=CMSw(i,j)+:.(i, j) 

8i829 


IFCrSMP,Cr,SECO'IO)30 to 86 

81839 


SEC 3 'IDs re 

91849 

86 

COMflNUS. 

91859 


SS(3-) = SeC3NO«C:NSWCiCONJ,J)-C(ICONJ, J) 

81860 

85 

CD^nWOE 

91879 


pMa=o, 

91889 


DO J38 

81890 


TE'iPsO, 

91999 


DO 8 3 iC = lfP 

91919 


IF(r£(£,K).LE.O)GO TO 88 

91929 


rE'<p-T£Hp^re(i,i^) 

91939 

88 


81949 


lFCrSMP.CiS.PMAX)G3 TD 188 

91959 


P»4A)tx,r£,»<p 

81969 



81979 

91989 

188 

continue. 

91999 

♦ ♦ 

♦ ♦ CHECKS THE FEASIBIUTY OF THE S063TI3M AMO JPOATSS THE 

92899 


INCJ'«8EMT IF NEGESSARY ♦♦ 

92019 



92029 

100 

lf(*'M0r.fSASE)G3 ro 17 

82039 


iCJsICJ^l 

82048 


iFCaCOSH-FIX^^A^J.Ge.IZOO) GO TO 19 

82058 


izoo«icasr+Fix\rAi4 

82069 


IfdCUGr.DGO T3 101 

92079 


CALGt RriMS;£S3) 

82089 


i.fe:ais^ii3-mi 

82099 


HRi:m23, 1110)1200 

•2tft9 

1110 

FORMiAffClXf 'THE FIRST FEASIBLE 50LUTI3N VALUE IS ',F8,2) 

i?ll9 


<IRmK23,llll)XFEAS 

P!i?t 

nil 

FORWClXr'TIME TAKEN FOR THE FIRST FEASIBLE SDLuriDN 

#1# 


1 HlEIi-LI. SECONDS') 

92149 


»RI.rS]U3,lil2)lCO0 

92159 

1112 

FORWfUr' NO, OF nodes EXAMINED FOR THE FIRST FEASIBLE 

92169 


1 SOLUTION ',16) 

92179 

lOl 

DO H J»1,N 

92189 

14 

IX(J.)=C0»EICJ,) 

92199 


ir(L6],GE:,£20D)G3 TO 16 , 

82299 

19 

KXCOOttO 

92219 


DO 130 Jal,N 

•222# 


lF(FLA3tJ))G0 T3 180 



gi224(> 


^xc<)(:3U)=j 

92259 

180 


92269 



92279 

♦ ♦ 

3R0S.RS THS. JOBS IN DECREASING ORDER 3F PENALTST J3LN0 HEAP SORT 

92289 


rFC<)cc,Du,6e,UG3 ra lei 

92299 


kk)c=!Cx::3u 

92390 


lFIRSrsKX:0U/2 

92319 

219 ' 

FlRSr^slFlRsr 

92329 

310 

JSFIRST 

92339 


ITE'^PSSXC J);rEMP5SE(rrEMP) 

92349 

320 


92359 ’ 


IF(^.Gr.KXX)GD ro 130 

92369 


IFC(<,[4T,<KX).AN0.(SECKXCK + 1)).L,T.(SE.(XX(i^))) JXsXfl 

92379 


lFCSS;CXXC<.)),Gr.'TE'lPjGO TO 330 

92389 


KX(J.)*<XCX) 

92399 


•SKi 

II 

92499 


30 n 320 

92419 

330 

fCXCJ)-l.rE4p 

92429 


FlRSrsFIRSr-l 

92439 


lFCFIKSr.3.r.0)G3 ro 310 

92449 


bXE^paXXCU 

92459 


XX(1)SKX(<XX) 

92469 


xxc<'XX)*i*rE'<p 

92479 


KXX»XKX-l 

92489 


IF(XXX,C0.'UG3 TO 181 

•2499 


IFIRST»1 

92599 


30 TO 219 

92519 



92529 

f« 

fixes the: jobs in the decreasing ORDER 3F PENALTI 

•2539 

181 

00 18. KKsl,KXCOU 

•2549 


J?««XCK1C:) 

•2559 


IF(F6<AG-CJ3) go to 18 



. lOPSffOPi-l 



ICDtI-*ICDUti 

•2589 


F t» A.Gt( J.) -T RO B • ' 

02599 


ICONjJsCDNElCJf) 

92689 


IS(r'3P)*lCD«J 

02619 


JSIXOPJSJ. 

•2620 


IVAli«(TQP)».l 

0263» 


DO 71 K?!l,P 

92649 

71 

pB(lC0NJ.,K,)apBCICO»J,K3-A(lCONJ,JfK) 

92659 


EL E 'i.EN aSIiSMEN -C NEN ( I CON J , J ) 

•2669 


KL3=XLB-:::(i:CaNJf J)-CNEWCICONJ#J) 

92879 


FIX'fALaFlXVAli+CCICONJiJ) 



9268(11 

18 


02699 


icosrteo.o 

027»9 


30 ro 16 

02719 

02729 

♦ » 

FrXS5 THE J38S 3f THE AGENT SEDECrEO FDR 8RANCR1M3. 

02730 

17 

SXCOjsD 

02740 


DD 23 J=l,N 

02760 


IF(!t(INDl, J).E0.O)3O ro 23 

02760 


IF(FW3.(J))G3 t3 23 

027 70 


!C;KE3j = iCXC3U + l 

02780 


ICXC!<x:31I)5J 

02790 

23 

:DNri^us. 

02899 


KKX*<X::3U 

02810 


iFCKXX.tie.DGD ro 35 , j 

02829 


IFI«ST»iCX:oU/2 1 

02830 

119 

FIRSrtelFlRST 

02840 

210 

jxfiRsr 1 

02850 


ITE4;P=KX( J.3;TEMPsSECITEMP) ! 

02860 

220 

K^2»j; ■ 1 

92870 


lF(X..Gr.KKX)G3 ro 230 j 

92880 


lF(U:.tir.XKX).AND.(SEtXX(X+l)) ,br.CSECKX(K)l ) ) ) XsXf 1 

92890 


lF(5!:!(KX(X)),GT,'rE><P)Ga TO 230 

92990 


KXC J3=XX(<) 

92919 


J»K 

02920 


so T'D 220 

•2930 

230 

KXCJi)ai:rE'<P 

92949 


prRSrteFXRST-l 

•2950 


lF(FIflsr.3.T.0)G3 ID 210 

02960 


t,TE:<P»KX(U 

02970 


KX(1)*KX(KKX) 

•2980 


XX(XKX3stilE«P 

92990 


KKx»xK:x-a 1 

* i 

i80i« 


IF(XXX,€Q.1)G3 ro 35 | 

•3019 


iFiRSTal' 1 



SO ro 119 

03039 

35 

iST’AiRs.lNoX: ■ 

03040 


00 28 Ji»l,KXCOU 

03060 


, rF(FOiA;S.CKXCX)5)SO rO 28 

03060 


DO 29 KssifP 

03070 

29 

lFCA;Ci;f«OI,KX(J),K).Gf.PS(lK01#K))GO TO 28 

03080 


JSfARSKXCJ;) 

03090 


roparop+t ' ■ ' , ' 

•3109 


ic.3Xs.rrou^i 

93110 


ISC np)ai STAR i 


831 3» 


lVALi(r3P)=l 

83148 


KCISTAR, jSrRR)=!l 

83158 


:0«SiCJSrAR)slST4R 

83158 


FLA3.C JSTAR) = .rRUE. 

83178 


DO 2 4 !C = 1,P 

85188 

24 

# 

PBC I.STiR,<)spBClSrAR,K)-A(lSTAR,jSrAR,K) 

83198 


fI)C\?AL = F.lXyALtCClSrAR,jSTAp) 

83298 


i(ti3=.a'3-CClSrAR, jSrAR)-CNEw(ISTAR,jSrAR) 

93218 


EtE'42NsEL2ME.'J-C'*E^(lSTAR,JSTAR) 

83228 


Ecasr-sizosT-ccisrAR, jSTAR) 

83238 

28 

:0Mri:«uE • 

83248 



83258 


:H5:iiCS THE J33S TO BE REAI.LUCATEO 

83258 

76 

iTRftCKa.FALSE. 

83278 


BO 7 7 f?(C = l,KXCO'-I 

83288 


JsK)t(K!C) 

83298 


IF (FtiAG.CJ)) GO T3 7 7 

83388 


• PMI'IaLARGS: 

83318 


SEC:3'lI5 = ti'ARGE 

83328 


BO 43 1 = 1, '4 

•3338 


iFCI.eO.lMODGO ro 48 

83348 


IF CFL(I,J)) GO ra 48 

83350 


DO 41. k:=i,,p 

83358 

41 

IF CMCI:,J,K).GT.PB(I,K)) GO TO 48 

83378 


rEMPxC'CI, J.)fCME<<CI,J) 

93360 


IF(rS)MP,GE;.P«lN) G3 to 43 

8339# 


SECOSOssPMifi 

83488 


pmin«.te:'<p . 

83418 


IKOCSEsi; 

83428 


GO ro- 48 

83438 

43 

If CfEMP.Gi:, SECOND) GO TO 48 

83448 


SEOOfO^XEW . 

83458 

48 

cdmumoe; 

•3458 


If (p.4!i;!i,c£;.0AftGE;) itrackb.true. 

83478 


XClUOCfJijaO 

83468 


X(IS0€X,J?*1 

83490 


SE-UiJaSEGOMD-PMlN 

83598 


:D»Si(j;)»I!40EX 

83518 

■, 

lE.OSm^DST-nClMDI,J)+CCINDEX,J) 

•3528 


BO 44 i<;sl,P 

83538 


tec i:4BEX,X;)«TEC INDEX, K)'>-AC INDEX, JrK) 

•3548 

44 

rE(i:'IBI,K)=TEClNDI,K)-AClNDl,J,IC) 

•3558 


rE« P =c;c I M BI , J 3 •♦•CNEW ( INDI , J> 

Jl35jl8 





gi357i) 

77 

:dmJ‘£mus 

»3539 

9359^ 


IF( I r'R<v:K)33 TO 16 

93699 

** 

:he:k' for pruning and selects .the agsnt fdr 8RR«:-iiNiG 

93610 

93 

FE&SSP.TRLJE. 

93620 


LB=<.L6-;SLSMEH + FU\^AL 

93630 


LBstiSfO, 99999 

93640 


ifcCi'b.ge.izdd) go ro is 

03659 


pMft)C5-lOOOO.O 

93660 


00 51 

03670 


00 5l <aifp 

93680 


IFCr'!::CI,K).Gr.O)FE/lSE=,FALSE. 

03690 

51 

:09r£>ij£ 

93700 


IF(F£A5E.)G0 to loa 

93710 


00 U5i r=i,H 

93720 


rSHPaO, 

03730 


00 151 K=l,P 

03740 


lF(l’S:(IfK).Le,0)G3 TO 151 

§3750 


rEHPsJeHPtTECIfK) 

•3760 

151 

councils 

93770 


lF(rs.MP,LS.PHAX)G3 TO 1151 

03789 


PMAICsTSHP 

0379o 


IN0£»I 

■3800 

llSl 

CDHTINUE. 

03810 

93820 

, 

30 TO 100 

93830 

♦ ♦ 

back: track 

93840 

16 

lF(nP.S0,O)G0 TO 78 

93859 


KKsl./A.LiCr3P) 

93869 


IslSCTOP) 

93879 


JfajSCTOP) 

93889 


If (KX:.€a,23GO T3 UO 

93899 


jyAWT0P)»2 

mk' 


'rL'CJ;#3t3«>TROE#' 

mnw- 


rLt!G<3.3*vfALSE, 

03920 


FIX V AIi*mf AL-C C I , J ) 

03930 


icosmeofit^xi^j.) 

93949 


00 120' 

03950 

120 

P8(i:, IC)*l>eci»K)tA£I,J,K) 

03960 


ELEH!ENseLEiMES+CMEW£I,0) 

93970 


pmi::'i«l-airge; 

' 93980 


5EC:3.'iO«LA'RG€ 

93990 


DO 81 IIsl^M 

94000 


IFCFLICII , J.))GO TO 81 



||4i02l> 

82 

lF(4(U,J,K),C;T.Pti(Il,K))G0 TD 81 

049301 


rE'<p=c(ii, 

04940 


lFCrSMP.GS..P'1IN)G3 TD 83 

84950 


Se:'3'l9sP'!l,M 

04960 


?m^=rs'ip 

04970 


iNOEiCslI 

04OB0 


30 n 81 

04090 

83 

IFCTSI'IP.GS.SECOMDIGO TO 81. 

04100 


SE^OOsTE'lp 

04110 

81 

aOMTlStlS 

04120 


lFCP'1l!i.G£..LARGE) GO TO 16 

04130 


ir3'lJ]=30.ME,C J) 

04140 


XdCOSJ, J)=0 

04150 


xci'Joex, J)=l 

04160 


SEC jd*sscond-pmim 

04170 

84 

iCl.B»'<li8 + P'irfJ 

04180 


iCOST^iroST+CCl^DEX, J)-C(IC0NJ, J) 

04190 


00 198 Kal,p 

04200 


r£.Ci:D?iJ.,fC) = rE(lC0NJ,K)-AClCONJ,J,K) 

04210 

108 

tec i:40EX,(C)=TE( INDEX, K)+AC index, J»K) 

04220 


2DNS]( J)*INDEX 

04230 


30 52 jsi,N 

04240 


If (FL.AG.( J)) G3 TO 52 

04250 


IFCFUICI, J)) GO TO 52 

04260 


IFCDNSXJO.EQ.I) GO TO 52 

04270 


00 54 XslrP 

04280 

54 

If (A,(I,J;,K.),GT,PBCI,K)) GO TO 52 

•4290 


irO'lJJsCDNSiC J) 

«430i 


TEMP^CCI, JDfCNE'JCI, J) 

04310 


rEMPl«ClCi:DMJ,J)+C?JEW(ICONj,J) 

04320 


IF C TEMPI. Oe, TEMP) GO TO 52 

04330 


XCI-QNj:,J)aO 

04340 


Kl»8*KtB;-!rEMPl+TEMP 

iiiii: 


OO 55 K>k1,p 

Ifilt 


irC I;, K) PTEt I,K)tACl,J,K9 

04370 

55 

miCDflJ;, fC.3aTEClCDMJ,K3-A(lCONJ,J,K) 

04380 

t 


tra s n?iCD s i -c c i :dm j., j ) +c c i , j 3 

04390 



04400 


XCX» JD»1 

04410 

52 

CDICT’IHUC 

04420 


30 TO- 93 

04430 

130 

FLCI, Ds.fAliSE, 

04440 


■ TOP-STOP.-! 

»*450 





044&OI 

»447» 

94480 

94499 

94599 

•4519 

•4529 

94559 

94549 

94559 

94569 

94579 

94589 

94599 

94699 

94619 

94629 

94639 

94649 

94659 

94669 

94679 

•4689 

94699 

94799 

•4719 

•4729 

94739 

94749 

94759 

94769 

94779 

94789 

i,47i9. 

■'■mm 

94819 

04829 

94830 

94849 

94859 

94860 
94870 
94880 
94890 


rS'^P=CMS,w(rC3f!j, j)fC(icoNJ,j) 
r£'<Pi=:;(i, j)4-CMSw(i,j) 
iFClTR^CiOGD ro 42 
lFCrS><Pl,3.T,TEMP)G3 xO 16 
42 XCICDHJ, j)=0 

Kt43=^L3-rSMP + rE'lPl 
£G3Sr«t3QSr + CC If JX(IC0NJ, J) 

SECJjsrS.MP-TEMpl 
30 124 KslfP 

rECtOD'iJfOsXEt rC3MJ,K)-A( iCONJfJfK) 

124 rECI,K)5.TS(I,K)VA(i, j,K) 

:0NS,( J )B1 
30 ro IS 

78 RriMEC'JS) 

irorALs'is-^Ni 
-#RirS!C23,102)rx0rAL 

102 P0R4AlTCU,'r3TAG TIME MlO) 

<iRirs:(23,l03)lC0U 

103 F0R4A!T(lXf* fDTAU MD OF NODES 'rl7) 

RErJRMjENO 

SOBROUriNS; GENRATCM, N,P,AfCfB,ISEED) 

** XHIS SUBRDuriME GENERATES A MULri-RESOUKCE GENSRAGiI SEO 

♦♦ ASSIGNMENT PRDbbE^ 

♦ ♦ if^**¥* **¥**¥ *¥*******$**** ********¥****•¥**¥$¥¥¥¥¥ ¥*¥***•**¥ 

dimension AC25,20I),6),CC25,200),B(25,6),S0M(25,5) 
iNrEGER P 

CALi.' SETRANCISESO) 

DO I I 3 1,4 
DO I K.s.lfP 
1 sumU;,k3so 

00 10 j;».ifN 
pHi:N«i 00000.0 

DO 2D IlslfM 
DO ID. K!s.i , p 
|K#j0#RAiM(X)t5 
C IP.PW:4^,99 

A'CIfO^KO^PK 
30 CON tote: 

PKsAO^rahCXJ+IO 

IPK'SPK. 

ClCI'f JDsIPK: 

lFC:!(I.,Ji?*<3^,PMlN) GO TO 20 



ei49ta 

20 

co^mus. 

•4929 


DO 4D K's.lfP 

•4939 

40 

SU'< c £?l DEX r Kl *SU«;C ttel^ r KHrIYnDEX , J , K ) 

•4949 

10 

:0Mrx:Mas, - ' ■; f'’ 

•4959 



•4969 


DO 53 1=1, M 

•4979 


DO 53 K=l,P 

•4989 


lFCSj'<(r,K,).LT,P'iA}C)GO TO 50 

•4999 


pMftKsSLJ4(I,K) 

85099 

50 

zomfimos 

•5019 


P8s0.6*M/'<1‘15<-0.4»PMAX 

•5028 


O 

ft 

95039 


DO 77 Kai,P 

•5049 

77 

BdrOaPB 

85059 


RETJRNjEND 



